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Abstract 


A method of dynamic substructuring is presented which utilizes a set of static Ritz vectors as a replacement 
for normal eigenvectors in component mode synthesis. This set of Ritz vectors is generated in a recurrence 
relationship, which has the form of a block-Krylov subspace. The initial seed to the recurrence algorithm is 
based on the boundary flexibility vectors of the component. This algorithm is not load-dependent, is 
applicable to both fixed and free-interface boundary components, and results in a general component model 
appropriate for any type of dynamic analysis. This methodology has been implemented in the 
MSC/NASTRAN normal modes solution sequence using DMAP. The accuracy is found to be comparable to 
that of component synthesis based upon normal modes. The block-Krylov recurrence algorithm is a series of 
static solutions and so requires significantly less computation than solving the normal eigenspace problem. 
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Introduction 

Component mode synthesis is a methodology for analyzing large structures by separating them into smaller 
components, which can then be recombined to analyze the entire system. This methodology has become well 
established and widely used in structural dynamic analysis. The advantages of component mode synthesis 
include, the lower computation costs associated with the smaller components which are analyzed, and the 
flexibility of data management gained by working with the discrete components. 

The typical component mode synthesis algorithm is briefly described [1]. A large structure is broken into 
components, with each component having a set of boundary, or interface, points. At these interface points, 
fixed or free boundary conditions are assumed, and a corresponding set of component normal mode shapes, 
or eigenvectors, is determined. The eigenvectors are augmented by a set of modes which are associated with 
the component’s boundary flexibility. Depending on whether a fixed or free interface is selected, these modes 
are the constraint modes or the attachment modes, respectively. The combined set of component normal 
modes and boundary modes are used to represent the component in subsequent system analysis, by using the 
following transformation process. The combined set of modes form a coordinate transformation matrix which 
transforms the physical coordinates of the structural model into a combination of modal coordinates and 
boundary coordinates. The boundary coordinates are retained in the physical space, so they can be used to 
couple the components for subsequent system analysis. 

A component’s size, although smaller than that of the complete structural model, can still be large enough 
to make computation expensive. The rapid reduction in cost per calculation in today’s digital computers has 
not necessarily led to a reduction in total computation cost. Instead, engineers have exploited the increased 
computational resources by creating larger structural and component models. The larger models have allowed 
for more structural details to be included, as well as more refined data recovery, but they may be expensive 
to formulate and analyze. In order to reduce the computational cost associated with large component models, 
it is desirable to develop more efficient methods of formulation. Since the solution of the normal eigensystem 
problem requires the largest computational effort in component formulation, it is logical to develop alternate 
methods which circumvent the eigensystem solution entirely. 

A method, which does circumvent the eigensystem solution, has been defined in literature and is briefly 
described [2-4] . The boundary flexibility modes, specifically either the same constraint modes or attachment 
modes that were mentioned previously, are multiplied by the component mass matrix to create a force matrix. 
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Static analysis is then performed, using this force matrix and the component stiffness matrix, to obtain a 
matrix, or block, of vector displacements. A recurrence relationship of matrix multiplications, which have been 
shown to be a Krylov sequence [3,5], then defines a series of matrices, or blocks, of vector displacements. The 
calculated vectors are orthogonalized, using normalized Gram-Schmidt orthogonalization [6], These vectors, 
which can also be thought of as static modes, replace the normal modes in the component formulation 
methodology. Because the vectors are calculated in blocks and are based on a Krylov sequence, the subspace 
defined by these vectors is called a block-Krylov subspace. 

The primary contribution described in this paper is the incorporation of the above methodology into 
MSC/NASTRAN. Since the block-Krylov vectors replace the normal modes in typical component mode 
synthesis, the physical to modal coordinate transformation, the constraint modes, and the attachment modes, 
are not changed. Therefore, the component mode synthesis formulation currently available in 
MSC/NASTRAN was utilized, except for the substitution of the normal modes by the block-Krylov vectors. 
In addition to the contribution described above, the technique of using Gram-Schmidt orthogonalization 
within each vector block was added to the methodology. 


Problem Definition and Theory 


Overview 

Wilson, Yuan, and Dickens [7] originally proposed the use of Ritz vectors, based upon external loading, for 
structural dynamic analysis. This formulation reduced an entire structure, not a component. The algorithm 
begins with a set of externally applied loads. The displacements from the static solution to the applied loads 
become the initial Ritz vector. That vector is then multiplied by the mass matrix to become the next force 
vector. This sequence is repeated to form a recurrence relationship. This recurrence relationship is used in the 
papers discussed below and throughout this work. 

Nour-Omid and Clough [5] investigated Wilson, et al.’s methodology and found that the proposed recurrence 
relationship was actually a Krylov sequence. A Krylov subspace of order j is a vector space defined by 

[^A^A 2 *,...,^" 1 *] to 


where <j> is a column vector and A is a square matrix. If A is n X n dimensional, and if j — n, the Krylov 
vectors span the n-dimensional space [3], and an exact solution can be produced. In structural dynamics, the 
Krylov subspace can be defined by the following sequence, 

[r,K Mr, (K -1 M ) 2 r, ... ,(K _1 M )* -1 r ] ( 2 ) 


where K is the stiffness matrix, M is the mass matrix, and r is a starting vector (or in block-Krylov, a set of 
vectors). This Krylov sequence is identical to the Lanczos eigenvalue extraction algorithm, when applied with 
complete orthogonalization with respect to both the mass and the stiffness matrices. 

The use of Krylov vectors was shown to be applicable to component mode synthesis by Wilson and Bayo [8]. 
The Ritz vectors calculated were based, once again, upon an external load. Only a formulation for components 
with fixed interface boundary conditions was presented. 

Craig and Hale [S], and Abdallah and Hucklebridge [2], demonstrated a methodology applicable to 
components with fixed or free interfaces, with or without rigid body modes, and with no applied loading. 
Components having no applied external loading were formed using the boundary flexibility matrix, multiplied 
by the mass matrix, to form a force matrix. This force matrix produces a set of vectors, which are referred 
to as a block. As discussed previously, the boundary flexibility matrix is defined as either the constraint modes 


3 


or the attachment modes, depending on whether Fixed or free interface conditions are selected. The 
methodology contained in these two papers is reviewed in the next three parts of this section. Abdallah and 
Hucklebridge also quantified the advantages, in computational effort, that block-Krylov vectors have over 
normal eigenvectors. 

Yiu and Landress [4] also developed a methodology for forming a component which does not have an external 
applied load. However, their formulation is applicable to fixed interface components only. A criteria for 
concluding the recurrence sequence, based upon the rigid body mass and flexibility represented by the 
calculated Eitz vectors, was proposed. 


Fixed Interface Methodology 


First, as is standard in component mode synthesis methods, the finite element component mass, m, and 
stiffness, k, matrices are partitioned into internal and external degrees of freedom, denoted by subscripts i and 
c, respectively. 


m cc 

m ci 

m ic 

m ii 


( 3 ) 


k = 




( 4 ) 


The constraint modes are defined by 


*!<■ -*«''** 


( 5 ) 


which is the same definition used in standard component mode synthesis. 

For Wilson’s method [7], a set of externally applied loads is required to obtain the initial set of Ritz vectors. 
For the boundary flexibility method, this set of loads is created by multiplying the constraint modes by the 
mass matrix. (Craig [3] also included the off-diagonal mass matrix in his formulation.) Since the mass matrix 
is used to create the loads, they can be considered inertia loads. This set of inertia loads are then used to 
generate the initial set, or block, of Ritz vectors using the following 


qi’= V^ii^ic^ic) 

where the superscript ** indicates that the vectors in the matrix have not been normalized. The first block 
of vectors is normalized using the following equation. The subscript r, in the following equation, signifies that 
the block is normalized vector by vector, and there are c vectors within each matrix, or block. 








r = 1,2,.. .,c 


( 7 ) 
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The subsequent sets of block-Krylov vectors are generated using the recurrence relationship [5,7], which was ' 
defined in Eq. (2), 


where the superscript * signifies that the vectors have not been orthogonalized or normalized. The additional 
sets of vectors are orthogonalized, with respect to the mass matrix, with all previous vectors. The process used 
to perform this orthogonalization is a normalized Gram-Schmidt procedure. 


where 

c=qi J _i T m ii q’ 

and is the concatenation of the previous sets of block-Krylov vectors, 


( 9 ) 

( 10 ) 


( 11 ) 


where &11 vectors have been normalized as follows. 


Vr 


‘‘J, 


/( V T -*0 


r = 1,2,..., c 


( 12 ) 


The complete set of calculated block-Krylov Ritz vectors is included in the transformation matrix as Q,. (The 
resulting transformation matrix has the same form as that of ” Craig-Bampton” component mode synthesis 
[9], with the Krylov vectors replacing the normal modes.) 


* = 




0 

Qi 


The physical mass and stiffness matrices are transformed into the component modal matrices 


li =1' T m'J' 
k =ty T k'I' 


(13) 


(14) 

(15) 


The resulting mass submatrices are 


where 


V- = 


^cc 

Mlc 


Md 


Mcc ■ $ ic T ( m ii $ ic + m ic ) +m ci*k + m c 


Mlc“M c i T =Qf(mi^ic + mic) 


(16) 


(17) 

(18) 
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(19) 


= % = ^1 ™uQi 


The resulting stiffness submatrices are 


where 


k = 


*cc *el 
*lc *11 


K cc ~ ^d^ic + ^cc 


*lc = *cl T = 0 

«n kaQi 


( 20 ) 


( 21 ) 

( 22 ) 

(23) 


The use of constraint modes in the transformation matrix leads to the null off-diagonal partitions of the 
component stiffness matrix, just as in the approach based upon normal modes. 


Free Interface Methodology for Components with No Rigid Body Modes 

When allowing the interface points of a component to be free to deflect while forming the component, a 
somewhat different basis for the initial vector of the Krylov algorithm is required. The attachment modes, 
rather than the constraint modes are utilized in initial block definition. By definition, the attachment modes 
are the columns of the flexibility matrix which correspond to the interface degrees of freedom. 


g^k’ 1 


(24) 


Ba 


gcc 

gic 


The initial block of vectors in the free interface formulation is defined as 




* * 


k _1 mg a 


(25) 


(26) 


and is normalized as follows. 


q x 

q, - i t - 1,2, (27) 

/(I,/ • T <°q,,") 

Note that the unpartitioned physical mass and stiffness matrices of the component are used in the free 
interface formulation. The recurrence algorithm then proceeds in the same manner as in the fixed interface 
methodology. 
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(28) 


V -k^mqi.j 


(29) 


c 


q^.l mqj 


(SO) 








r = 1,2,..., c 


(31) 


Formation of block-Krylov component then follows the normal component mode synthesis techniques which 
were presented by MacNeal [10] and Rubin [11]. To combine the ”Rubin-MacNeal* method with the presented 
method, the normal eigenvectors are simply replaced with the block-Krylov Ritz vectors, as in the fixed 
interface methodology. The free interface methodology uses residual flexibility terms, which fully define the 
stiffness missing from the modal space due to excluded modes and are described below. The flexibility 
contained in the calculated Krylov vectors is given by the following equation. 

8 k = Qi(Qi Tk Qi)“ 1 Qi T (32) 


The unrepresented flexibility, or residual flexibility, is defined as 


8<3 ” % ~ 8k 


(SS) 


The residual flexibility matrix is then partitioned in the same manner as the flexibility matrix was in 
Eq. (25), when the attachment modes were created for initial vector calculation* The result is the residual 
attachment modes* 




a 



(34) 


When the residual attachment modes, g a<J , are added to the Krylov modes, Qj, the complete flexibility of the 
component is represented. 

The residual attachment modes and the Krylov modes are used to form the component transformation matrix. 
This matrix transforms the physical subspace, u, to the modal subspace, p, and is defined by the following 
equation. 




* 

£> 

Pc 

V 


8ic d Oil 

.PI. 


(35) 


In order to provide physical interface degrees of freedom, for use in component coupling, p c in the above 
equation is back-transformed to eliminate it from the right-hand side of the equation. This results in the 
following transformation matrix, 
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( 36 ) 


where 


«c 


Icc 

0 



.V 


3ic 

Qa*. 


.Pi. 


®>ic Kjc^Kcc^ 


-1 


(37) 


Q/ =Qil-8ic d 8cc d _1 Qcl 


(38) 


The transformation of the component mass and stiffness matrices then proceeds in a similar manner as shown 
in Eqs. (14) to (23), with the following differences. The Q a matrix partition replaces the Qj matrix partition. 
The g ic * matrix partition replaces the $ ic matrix partition. In the fixed-interface methodology, the definition 
of the constraint modes, $ ic , leads to terms in the component stiffness matrix which cancel out. In the free- 
interface methodology, the definition of the transformation submatrices has changed and so this cancellation 
does not occur. Therefore, Eqs. (21) and (22) are replaced by the following equations. 


and 


*cc 52 8ic (^iiKic *^ic) + ^ci8ic + ^ c 
*jc = *cj T ■ Qa* T (Mic* +k ic) 


(39) 

(40) 


Free Interface Methodology for Components with Rigid Body Modes 

When a component has rigid body modes, the associated stiffness matrix is singular. The inverse of the 
stiffness matrix, the flexibility matrix, cannot be directly obtained, and therefore the attachment modes 
cannot be directly obtained. To circumvent this problem, Rubin [11] presented the following method for 
obtaining the residual elastic attachment modes of a component with rigid body modes. 

First, the stiffness matrix is constrained from rigid body motion by partitioning out r degrees of freedom, 
where r is the number of rigid body modes. 


k = 



kwr 


The remaining partition is then inverted. 


&WW = ^ww 


This flexibility matrix is then expanded back to n (w + r) size. 


®wr 



(41) 


(42) 


( 43 ) 


A square projection matrix is defined by 


8 


A-I„n-mi r t r T 


(44) 


where $ r is the rigid body modes matrix. The elastic flexibility matrix, g € , with rigid body motion removed, 
is shown in reference [ll] to be 


g e = A T g c A ( 45 ) 

Now the analysis proceeds in a similar fashion to the previously discussed methodology of the free interface 
component with no rigid body motion. The major difference between the two approaches is that the elastic 
flexibility matrix is used in place of the general flexibility matrix. The inertia relief attachment modes are 





(46) 


The initial block of vectors is calculated using the inertia relief attachment modes and the elastic flexibility 
matrix. 


qi” 


(47) 


The subsequent block-Krylov Ritz vectors are calculated, orthogonalized, and normalized as shown in 
Eqs. (28) to (31). The residual elastic flexibility terms are also calculated as shown in the free interface with 
no rigid body modes discussion, Eqs. (32) to (34). Creation of the transformation matrix, Eqs. (35) to (40), 
is also similar to when no rigid body modes are present. The one exception is that the rigid body modes must 
be included in the transformation matrix. Therefore, Eq. (35) is replaced by 



Qcl $ cr 
Qil *ir 


Pc 

Pi 

Pr 


(48) 


Formation of the final transformation matrix, and subsequently the component mass and stiffness matrices, 
is then performed as described in the previous section. 


Discussion 


As previously discussed, most component mode synthesis applications use the normal eigenvalues of the 
substructure to form the component. It is accepted that if all the eigenvalues of a system are used to form the 
component, an * exact” fmite element solution may be obtained. The same fact holds true for components 
based upon block-Krylov subspaces. Mathematically, this is proven in references [3] and [12]. If n orthogonal 
block-Krylov Ritz vectors are used to form a component of a n-size system, the same "exact” solution as from 
normal eigenvectors is obtained. 


Orthogonality 

Upon completion of the block-Krylov vector calculations and Gram-Schmidt orthogonalizations, all vectors 
are re-checked for orthogonality. Orthogonality checking is required because the normalized Gram-Schmidt 
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orthogonalization algorithm is not always successful at producing independent vectors [6], The Gram-Schmidt 
procedure will fail on occasions when vectors, although theoretically independent, are dependent within the 
numerical constraints of current digital computers. These vectors cannot then be made orthogonal and 
independent, using the normalized Gram-Schmidt procedure described in the previous sections. 

The orthogonality of the block-Krylov vectors is re-checked using the following 

L = Qi T m ii Qi ( 49 ) 


If all vectors are orthogonal and normalized with respect to the mass matrix, L will be a 1-size identity matrix. 
The mass matrix used in Eq. (49) is appropriate for fixed interface modes. For the free interface approach, 
the full physical mass matrix is used. 

Several solutions to the numerical dependence problem are being investigated. These potential solutions 
include selective re-orthogonalization, the modified Gram-Schmidt procedure [6], and Kahan’s 
orthogonalization procedure [IS]. 

During examination of the boundary flexibility method, it was found that the Ritz vectors within the Krylov 
blocks were dependent on each other. In fact, it was determined that there is no theoretical basis why the 
vectors within the initial block should be independent. The boundary flexibility algorithm, as presented in 
reference [2] and discussed in section two, makes no orthgonality check of the vectors within the Krylov block. 
Furthermore, the subsequent Gram-Schmidt procedure is ineffective because the blocks are orthogonalized 
with respect to a set of vectors that are not orthogonal. 

The problem of dependent vectors within the block can be eliminated by modifying the boundary flexibility 
algorithm. In order to make all the vectors in the solution orthogonal, each Krylov block was partitioned into 
it’s individual vectors. The Gram-Schmidt orthonormalizing procedure, shown in Eqs. (9) to (12), was then 
applied vector by vector to each Krylov block. The constituent vectors are then reassembled and the analysis 
proceeds as before. An alternate solution to the problem of dependent vectors within the blocks, is shown in 
reference [12]. In this solution a singular value decomposition is performed on the q- T mqj subspace. The 
resulting similarity transformation orthogonalizes the block. 


Modal Selection 

A subject which requires future investigation is the selection of block-Krylov component modes. The Krylov 
vectors are modes, although they are static modes (as are constraint and attachment modes), in contrast to 
the normal modes based upon the eigenvalues. The simplest approach to modal selection is modal truncation, 
based upon the numerical value of the natural frequency. However, modal truncation cannot be used directly 
with Krylov modes, because there are no eigenvalues associated with them. Modal selection techniques more 
sophisticated than modal truncation will also be investigated. The algorithm currently implemented has no 
modal selection capability. The Krylov modes produced are the Krylov modes used, in component formation. 


Example 


Programming 


The previously described algorithms were implemented in MSC/NASTRAN [14]. The use of a standard, 
commercially available computer program allows the results of this work to be easily used by other structural 
dynamists. Adding these methodologies to MSC/NASTRAN is allowed through the use of the internal 
programming language called DMAP (Direct Matrix Abstraction Programming). The standard solution 
sequences of MSC/NASTRAN are written in DMAP, and the source code of MSC/NASTRAN is available 
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at the DMAP level. For example, Eqs. (3) to (5), and Eqs. (16) to (23) are currently contained in the 
standard MSC/NASTRAN normal modes solution. Equations (6) to (12) were coded using DMAP, and were 
then incorporated into the MSC/NASTRAN solution sequence. A DMAP listing of the present methodology 
is contained in the Appendix. 


Definition 

An example case of a cantilevered beam was derived from a finite element model of the Space Station 
Freedom photovoltaic array central mast. The length of the beam was 1179.9 in. The modulus of elasticity, 
E, was 10.1 Xl0 6 lb/in. 2 and the moment of inertia of the cross section was 108.9 in. 4 Its weight per unit 
length was 0.2296 lb/in. The cantilevered beam was modeled with eleven nodes and 10 beam finite elements. 
A variety of boundary conditions, described in the next section, were imposed upon this beam. 


Results 

A component representation of the 10 element beam was created using the boundary flexibility method with 
block-Krylov iteration. The fixed interface approach, with two Krylov blocks and constraint modes, was used 
to form the component. The interface of the component consisted of one node and 6 degrees of freedom. The 
number of constraint nodes is equal to the number of interface degrees of freedom, and the size of the Krylov 
block is equal to the number of constraint modes. Therefore, each Krylov block contained six vectors. Since 
the component was formed with two Krylov blocks, it contained a total of 12 generalized coordinates. 

Plots of the lateral Krylov vectors, which represented the cantilevered beam, are shown in Figs. 1 to 4. The 
unorthogonalized vectors, as output by Eqs. (6) and (8) are shown in Figs. 1 to 2(a) and Figs. 3 to 4(a), 
respectively. The first normalized vector, as output by Eq. (7), is given in Fig. 1(b). (The first vector does 
not need to be orthogonalized.) The remaining orthogonalized and normalized vectors, as output by Eq. (12), 
are given in Figs. 2 to 4(b). The first two unorthogonalized Ritz vectors (plotted in Figs. 1 to 2(a)), which 
are in the first Krylov block, appear to be nearly identical. The first mode is similar to the classic first 
bending normal mode shape of a beam and, after Gram-Schmidt orthogonalization, the second vector has 
become the classic second bending normal mode shape (shown in Figs. 1 to 2(b)). 

The eigenvalues of the reduced Krylov subspace were calculated next. The first five natural frequencies from 
this reduced system are shown in Table 1. For comparison, Table 1 also includes the frequencies of a reduced 
system where the component was formed using traditional normal modes. This component was also formed 
with a fixed interface, but thirteen normal modes were used for numerical convenience. The full, or ” exact,” 
finite element eigenvalue solution is also shown. In the case of the Krylov vectors, no modal selection of any 
kind was used. For the case of the normal modal component, modal selection by truncation was used. The 
superior accuracy of the normal modal component, in the fourth bending mode, does not necessarily represent 
a limitation of the block-Krylov method, but instead demonstrates the need for a Krylov modal selection 
criteria. 

In addition to the fixed interface example, two free interface examples were created. Both were was based on 
the same 10 element beam described above, but with different boundary conditions. The first model was free- 
fixed, with the component interface being the free boundary condition, and hence it had no rigid-body modes. 
Equations (24) to (40) define the formulation of this free interface component. Table 2 shows the first five 
system natural frequencies of this model compared to the natural frequencies of the full finite element model. 
The other model consisted of the 10 element beam with free-free boundary conditions and rigid body modes. 
Equations (41) to (48) define the formulation of this component. Table 3 shows the first four elastic 
frequencies of this model compared to the frequencies of the full finite element model. In the two free interface 
cases, there is no comparison with the normal component mode synthesis. This is because standard free- 
interface MSC/NASTRAN routine does not use the "Rubin-MacNeal” method, and so a direct comparison 
was not performed. 
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Conclusions 


Using boundary flexibility matrices to initialize the block-Krylov recurrence algorithm provides an efficient 
and simple method for generating static Ritz vectors. Static Ritz vectors so generated accurately represent 
the dynamics of a substructure. Because this methodology does not require the solution of the component 
eigenvalue problem, the component can be formed with a significant decrease in computational effort. 
Although the issues of vector dependency and modal selection require further investigation, the block-Krylov 
Component Mode Synthesis method is a promising alternative in dynamic substructuring. 
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TABLE 1.— COMPARISON BETWEEN THE CANTILEVERED BEAM 
FREQUENCIES USING FIXED INTERFACE COMPONENT 
NORMAL MODES, KRYLOV MODES, 

AND FINITE ELEMENTS 


System frequencies 
full finite element 

From component 
normal modes, 
A% 

From component 
Krylov modes, 
A% 

1st Bend 

0.5464 Hz 

0.5464 Hi 

n 

0.5464 Hi 

wm 

2nd Bend 

3.415 Hz 

3.415 Hi 

H 

3.415 Hi 


3rd Bend 

9.522 Hz 

9.522 Hi 

H 

9.610 Hi 


4th Bend 

18.57 Hz 

18.57 Hz 

rf 

24.24 Hi 

mm 

1st Tors 

27.58 Hz 

27.58 Hz 

0 . 

27.58 Hi 

wtm 


TABLE 2.— COMPARISON BETWEEN THE CANTI- 
LEVERED BEAM FREQUENCIES USING FREE 
INTERFACE COMPONENT NORMAL 
MODES, KRYLOV MODES, AND 
FINITE ELEMENTS 


System frequencies 
full finite element 

From component 
Krylov modes, 
A% 

1st Bend 

0.5464 Hi 

0.5464 He 

0 . 

2nd Bend 

3.415 Hi 

3.415 Hi 

0 . 

3rd Bend 

9.522 Hz 

9.524 Hs 

.02 

4th Bend 

18.57 Hz 

19.23 He 

3.6 

1st Tors 

27.58 Hz 

27.S8 He 

0 . 


TABLE 3.— COMPARISON BETWEEN THE FREE- 
FREE BEAM FREQUENCIES USING FREE 
INTERFACE COMPONENT KRYLOV 
MODES, AND FINITE ELEMENTS 


System frequencies 
full finite element 

From component 
Krylov modes, 
A% 

Rigid Body (6) 

0. Hz 

0.000 Hz 

0 . 

1st Bend 

3.474 He 

3.474 Hz 

0 . 

2nd Bend 

9.549 Hz 

9.549 He 

EH 

3rd Bend 

18.65 Hi 

18.69 Hz 

Mm 

4th Bend 

30.70 Hi 

33.15 Hz 

ESI 
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(a) Initial calculation. 



(b) After normalization. 


Figure 1 . — Fixed interface Krylov mode, block one, vector 
one. 



(a) Initial calculation. 



(b) After Gram-Schmidt orthogonalization. 

Figure 2.— Fixed Interface Krylov mode, block one, vector 
two. 



(a) Initial calculation. 



Figure 3.— Fixed Interface Krylov mode, block two, vector 
one. 



(a) Initial calculation. 



(b) After Gram-Schmidt orthogonalization. 

Figure 4.— Fixed interface Krylov mode, block one, vector 
two. 



Appendix - MSC/NASTRAN DMAP Listing of the Boundary 
Flexibility Method of Block-Krylov Component Synthesis 


$ 

S BOUNDARY FLEXIBILITY METHOD DMAP IMPLEMENTATION 
$ 

9 VERSION 1.0 - KELLY CARNEY 

9 IMPLEMENTED ABDALLAH’S WORK AS PUBLISHED 

$ 

$ NOTES ON USAGE 
S 

$ - A DUMMY METHOD AND EIGR CARD MUST BE USED TO INITIATE RITZ 
$ VECTOR CALCULATION. 

S 

$ - MODAL DEGREES OF FREEDOM MUST BE PROVIDED AND SPECIFIED 
9 ON SEQSET CARDS, AS WITH STANDARD NASTRAN SOLUTION 

$ 

$ - MIXED BOUNDARY CONDITIONS (B AND C SET POINTS) ARE NOT YET 
S IMPLEMENTED. USE B OR C, NOT BOTH. 

9 

9 - NOT COMPATIBLE WITH DYNAMIC REDUCTION, IF REQUESTED IT WILL 
9 NOT BE PERFORMED 

S 

I - AUTO-OMIT FEATURE HAS BEEN DISABLED. IF AUTOSPC IS ON, 

9 AUTO-OMIT IS NOT REQUIRED 
$ 

$ PARAMETERS 

t 

$ PARAM,BFLEX,-1 - A NEGATIVE VALUE ON THIS PARAMETER TURNS ON THE 
9 BOUNDARY FLEXIBILITY METHOD. 

$ IF MULTIPLE SUPERELEMENTS ARE BEING PROCESSED, 

$ PLACE IN SUBCASE SECTION OF CASE CONTROL 

$ 

9 PARAM,L,I - L IS THE NUMBER OF RITZ VECTOR DESIRED. IT IS 
$ SET BY THE NUMERICAL VALUE OF I 

$ 

S PARAM,DIAG,1 - A POSITIVE VALUE CAUSES A LARGE AMOUNT OF MOSTLY 
9 1NDESCIPHERABLE DATA TO BE PRINTED 

$ 

$ MATRICES 

9 

9 INPUT - CMKXX, CMMXX, GOAT, USET 

% 

$ OUTPUT, FIXED INTERFACE - GOQ, 

$ OUTPUT, FREE INTERFACE - GOQ, REVISED MAA AND GOAT 
$ 

9 

9 SETUP BOUNDARY FLEXIBILITY PARAMETER AND DISABLE DYNAMIC 
9 REDUCTION 
9 

ALTER 720 I 

TYPE PARM,,I,Y,(BFLEX=I) S 
IF( BFLEX < 0 ) NODYNRED = -1 9 

S 

$ DISABLE AUTO-OMIT REDUCTION 

9 

ALTER 767 9 

IF( BFLEX < 0 ) NOARED = -1 $ 

9 

9 MOVE TO BOUNDARY FLEXIBILITY SECTION IF REQUESTED 

S 

ALTER 775 9 

COND BOUNFLX,BFLEX $ JUMP TO BOUNDRY FLEXIBILITY DMAP SECTION 
S 

ALTER 763 $ 

JUMP ENDBFt 
9 

9 BEGIN BOUNDARY FLEXIBILITY METHOD SECTION 
9 

LABEL BOUNFLX S 

I 

S CREATE AND SET INITAL PARAMETERS 

S 


TYPE PARM„I,N,(I=2) 9 
TYPE PARM,,I,N,(Il=2) 9 
TYPE PARM,,I,N,(J=2) S 
TYPE PARM,,I,N,(IR=2) $ 

TYPE PARM„I,N,(JSET) S 


INCREMENT OF RITZ VECTORS 
INC OF RITZ VEC IN FIRST SET 
INC OF RITZ VEC IN ITH SET 
INC OF R-SET 
SIZE OF RITZ VEC SET 
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TYPE P A RM„RS,N, (MAXRAT) 3 MAXRAT PROM DECOMP OF CMKXX 

TYPE PARM m I,Y,(L=1) S TOTAL NUMBER OF RITZ VEC 

TYPE PARM„l,Y i (DIAG=-l) * DIAGONOSTIC PRINTOUTS 
FILE QP=APPEND,OVRWRT $ ALLOW QP TO APPEND, OVERWRITE 

FILE QL=APPEND 3 ALLOW QL TO BE APPENDED TO 

FILE QNEW=APPEND S ALLOW QNEW TO BE APPENDED TO 

FILE PHlVZ=OVRWRT $ ALLOW PHIVZ TO BE OVERWRITTN 

FILE NPBIVR= APPEND 3 ALLOW NPHIVRTO APPEND 

FILE MAA=OVRWRT S ALLOW MAA TO BE OVERWRITTEN 

FILE KAA=OVRWRT $ ALLOW KAA TO BE OVERWRITTEN 

$ 

$ CHECK TO SEE IF BOTH B AND RC HAVE BEEN DEFINED 
S 

IF(NORC>0 AND NOBSET>0) MESSAGE / 

/’FATAL ERROR - BOTH BSET AND CSET OR RSET DEGREES OF FREEDOM’/ 
’HAVE BEEN DEFINED, ONLY ONE TYPE IS CURRENTLY ALLOWED’/ * 
IF(NORC>0 AND NOBSET>0) JUMP RFERR S 
S 

$ PARTITION OUT SUPORT DEGREES OF FREEDOM, IF REQUESTED 
$ 

VEC USET/VVRO/’V’/’R’/’COMPV $ CREATE PARTN VECTOR R-SET 

IF(NORSET >0) THEN $ EXECUTE ONLY IF R-SET 

PARTN CMKXX, VVRO,/„,K22/-l/////e 3 PARTION STIFFNESS MATRIX 

S 

ELSE $ IF NO R-SET 

EQUIV CMKXX, K22/ALWAYS S RENAME CMKXX TO KZZ 

ENDIF S 
3 

$ CALCULATE INITIAL RITZ VECTORS MATRIX 

S 

DECOMP KZZ/LZZ,/1/////S,N, SING/S, N.NBRCHG/ 

S,N, MAXRAT/ $ DECOMP STIFFNESS MAT 

IF (SING=-1 OR NBRCHG>0 OR MAXRAT> l.ES) MESSAGE/ 

/’FATAL ERROR-STIFFNESS SINGULAR-SUPPORT DOFS MAY BE REQUIRED’/S 
IP (31NG=-1 OR NBRCHO0 OR M AXR AT> 1 -E5) JUMP RFERR S SINGULAR MAT 
DIAGONAL KZZ/IZZ/’SQUARE’/O. $ Z-9IZE IDENTITY MAT 

FBS LZZjJZZ/KINV/ S CALCULATE K**-l 

3 

$ FOR FIXED BC METHOD 
3 

IF(NOBSET>0) THEN $ 

MPYAD KINV.CMMXX./KINVM/ $ CALCULATE K**-1*M 
MPYAD KINV,MOA,/KINVMOA/ $ CALCULATE K"*-l*MOA 
MPYAD KINVM,GOAT,KINVMOA/QlS/ $ CALCULATE INITIAL VEC 
UPARTN USETjQIS/QISA.m/’A’/’BV’Q’/J S PARTN DOWN TO B-SIZE 
ENDIF 3 

$ 

3 THIS SECTION IS EXECUTED FOR FREE INTERFACE METHOD 

S 

IF(NOROO) THEN 3 IF EITHER R OR C SET 

S 

$ THIS 9ECTION IS EXECUTED FOR FREE INTERFACE METHOD, 

$ WITH RIGID BODY MODES 

$ 

IF(NOR9ET>0) THEN 3 ONLY IF R SET 

$ 

S CREATE AND ORTHO-NORMALIZE RIGID BODY TRANSFORMATION MATRIX 

S 

VECPLOT, ,BGPDTS,EQEXINS,CSTMS„USET/PHIGR/GRDPNT//« $ 

UPARTN USET,PH1GR/PHIVR JI> /’G , /’VV’SVI 3 RIGID BODY MAT 

s 

$ EXTRACT FIRST VECTOR FROM MATRIX AND NORMALIZE 
3 

MATMOD PHIVR^n./PHIVRiyi/l $ EXTRACT 1ST COLUMN 
SMPYAD PHIVR1,CMMXX,PHIVR1, > ,/NPH1/S////1////1S 
DIAGONAL NPHl/NPHlS/*SQUAREV-.5 « SCALE FACTOR 
MPYAD PHIVRl,NPHlS,/NPHIVR S NORMALIZE VECTOR 
S 

3 ORTHOGONALIZE REMAINING RIGID BODY VECTORS 

3 

DO WHILE (IK <= NORSET) S 
3 

MATMOD PHIVR IJJ „/PHIIR,/1/IR * EXTRACT IRTH COLUMN 

I 

$ ORTHOGONALIZE IR VECTOR WITH PREVIOUS VECTORS 

3 

SMFYAD NPHIVR ) CMMXX,PHIlR ) „/CIR/3////l S SCALE FACTOR 
MPYAD NPHrVR,CIR,PHIIR/PHHRS//-l $ ORTHOGONALIZE 
3 

3 NORMALIZE ORTHOGONAL IR VECTOR 

3 
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SMPYAD PHIIRS,CMMXX,PHIIRS,,,/NPHIIR/3////l////l$ 

DIAGONAL NPHIIR/NPHIIRS/’SQUAREV-.S $ SCALE FACTOR 
MPYAD PHIIRS,NPHIIRS,/PH1IRSS $ NORMALIZE VECTOR 

$ 

APPEND PHIIRSS,/NPHIVR/2 s add new vector 

s 

IR = IR + 1 S INCREMENT 

ENDDO S END DO LOOP 

s 

S FORM INERTIA RELIEF MATRIX, GE (EQNS 21,22,23 - ABDALLAH) 

S 

MERGE, lls KINV,VVRO,/GC/-l//6 S EXPAND KINV TO X-SIZE 
DIAGONAL GC/IXX/’SQUARE’/O. S CREATE IDENTITY MAT 
SMPYAD CMMXX,NPHIVR,NPHIVR,IXX„lXX/AXX/4/-l/l////l/ * 

SMPYAD AXX,GC,AXX,„/GE/S////l////« $ CALCULATE GE 
S 

$ THIS SECTION IS EXECUTED IF FREE INTERFACE WITHOUT RIGID-BODY 
$ 

ELSE $ 

EQUIV KIN V,GE/AL WAYS $ IF NO R-SET 

ENDIF S 

$ 

$ FOR FREE INTERFACE, RIGID BODY OR NOT 
S 

VEC USET/WCO/’V'/’CV'COMP 1 / $ C- SET PA RTN MATRIX 
ADD VVRO,VVCO/VVRCO///l $ R AND C-SET PARTN MAT 
PARTN GE,VVRCO,/GVRC,„/l/ S PARTITION OUT O COLS 
MPYAD GE,CMMXX,/KINVM/ S CALCULATE K**-1*M 

MPYAD KINVM,GVRC,/QlSA/ $ CALCULATE INITIAL VEC 

ENDIF S 


$ EXTRACT FIRST VECTOR FROM MATRIX AND NORMALIZE FOR FIXED OR FREE 
S 

MATMOD QlSA, m ,/QlPS,/l/l $ EXTRACT 1ST COLUMN 
SMPYAD Q1PS > CMMXX,Q1PS > „/NQ1/S////1////1S 
DIAGONAL NQl/NQlS/'SQUAREV-.5 * SCALE FACTOR 
MPYAD Q1PS,NQ1S,/QP S NORMALIZE VECTOR 


S 

s 

COND PRINT1 ,DIAG S 

MATPRN KZZ,CMMXX,GOAT,PHIGR,PHIVR/ $ 
MATPRN NPHIVR,GC,AXX,GE,KINVM/ $ 
MATPRN GVRC,Q1SA,Q1PS,,/ 9 
PRTPARM ///"NOBSET'/l S 
PRTPARM ///'NOTSET’/l S 
PRTPARM ///’NORCVl $ 

PRTPARM ///' L’/l f 
LABEL PRINT1 S 
S 

S ORTHOGONALIZE FIRST SET OF VECTORS 
$ 


DO WHILE (11 <= NOTSET) $ 

S 

MATMOD QlSA,,,,,/Qlll ,/l/Il S EXTRACT HTH COLUMN 

S 

% ORTHOGONALIZE II VECTOR WITH PREVIOUS VECTORS 

S 

SMPYAD QP,CMMXX,QlIl,,,/CIl/3////l S ORTHO SCALE FACTOR 
MPYAD QP,CU,QlH/QmS//-l $ ORTHOGONALIZE 

S 

S NORMALIZE ORTHOGONAL 11 VECTOR 

S 

SMPYAD QlIlS,CMMXX,QlIlS,„/NQlIl/3////l////lf 
DIAGONAL NQlIl/NQlIlS/’SQU ARE*/-.5 S SCALE FACTOR 
MPYAD Q1HS,NQ1I1S,/Q1I1SS $ NORMALIZE VECTOR 

$ 

COND PRINT2.DIAG S 
PRTPARM ///» 

MATPRN QP,QlIl,CIl,QlIlS,qiUSS/ $ 

LABEL PRINT2 < 

S 

APPEND QlUSS,/QP/2 S ADD NEW VECTOR 

S NOTE: IF NEW VECTOR IS NULL IT IS NOT APPENDED 

S 

II = II + 1 S INCREMENT 

ENDDO S END DO LOOP 


$ 

$ INITIALIZE QL MATRIX 
S 

COPY QP/QL/ALWAYS/1 S QL - SUM OF VECTORS 

$ 
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COND PRINTS, D1AG S 

SMPYAD QL l CMMXX,QL,„/OCHKll/3////l S 
MATPRN QL.OCHKll,,,/ S 
LABEL PRINTS $ 

S 

$ BEGIN DO LOOP OF SOLVING FOR SUBSEQUENT VECTORS 

% 

DO WHILE (I <= L) $ 

$ 

$ CALCULATION OP NEW SET OF VECTORS 
$ 

MPYAD KINVM,QP,/QIS/ $ (K**-1)*RITZ PREVIOUS 

PARAML QIS/ /’TRAILER’/l/SjNjNOJSET $ NO. COLUMS IN QIS 

S 

$ MAKE NEW SET OF VECTOR ORTHOGONAL TO PREVIOUS SETS 

S 

SMPYAD QL,CMMXX I QIS„,/CI/S////1 I CALCULATE SCALE 
MPYAD QL,CI,QIS/QISS//-1 t ORTHOGONALIZE 

f 

$ EXTRACT FIRST VECTOR FROM NEW SET AND NORMALIZE 
S 

MATMOD QISS,,,,,/QIJl,/l/l $ EXTRACT 1ST COLUMN 

SMPYAD QIJl,CMMXX,QIJl > „/NQU/3////l////l$ 

DIAGONAL NQI1/NQI1S/‘SQUARE7-.S $ SCALE FACTOR 
MPYAD QIJl,NQIlS,/QNEW I NORMALIZE VECTOR 

$ 

COND PRINT4.DIAG $ 

PRTPARM ///T/l/S 
MATPRN QIS,CI,QISS,QIJ1,QNEW/ S 
LABEL PRINT4S 
S 

J = 2 % RESET LOOP COUNTER 

$ 

t ORTHOGONALIZE NEW SET OF VECTORS WITH EACH OTHER 

S 

DO WHILE (J <= NOJSBT) 3 

S 

MATMOD QISS„,„/QlJ,/l/J * EXTRACT JTH COLUMN 

* 

$ ORTHOGONALIZE JTH VECTOR WITH PREVIOUS VECTORS 

» 

SMPYAD QNEW ,CMMXX ,QI J,„/C J/3////1 S ORTHO SCALING 
MPYAD QNEW,CJ,QIJ/QIJS//-1 S ORTHOGONALIZE 

$ 

C NORMALIZE ORTHOGONAL JTH VECTOR 

$ 

SMPYAD QIJS,CMMXX,QUS m ,/NQIJ/S////1////1S 
DIAGONAL NQIJ/NQIJS/*SQUARE7-.S S SCALE FACTOR 
MPYAD Q1JS,NQ1JS,/QIJSS t NORMALIZE VECTOR 

$ 

APPEND QIJS3./QNEW/2 $ ADD NEW VECTOR 

S 

COND PRINTS, DIAG $ 

PRTPARM ///‘jyiS 

MATPRN QIJ.CJ.QIJS.QIJSS.QNEW/ S 
LABEL PRINTS $ 

S 

J = J + l $ INCREMENT 

ENDDO S END DO LOOP 

S 

t APPEND NEW VECTOR TO VECTOR MATRIX AND COMPLETE LOOP 

S 

APPEND QNEW./QL/2 S APPEND NEW VECTOR TO PREVIOUS 
COPY QNEW/QP/ALW AYS/1 S RESET PREVIOUS VECTOR TO NEW 

S 

I = I + 1 s INCREMENT 

ENDDO S END DO LOOP 

$ 

% FINISH PROCESSING MODE SHAPE MATRIX 

S 

IF(NORSET>0) THEN $ IF R-9ET, APPEND 

APPEND NPHIVR,QL/CMPHIXZ/1 $ RIGID BODY MODES 
ELSE $ IN NO R-SBT, JUST 

EQUIV QLjCMPHIXZ/ALWAYS $ CONVERT TO NASTRAN NAME 

ENDIF * 

S 

$ PERFORM MODE ORTHOGONALITY CHECK 

S 

SMPYAD CMPHIXZ,CMMXX,CMPHIXZ jn /OCHKLL/3/ ///l S PHI-T • M • PHI 
SMPYAD CMPHIXZ,CMKXX,CMPHIXZ,„/KBAR/S////1 S PHI-T • K • PHI 

$ 
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DIAGONAL OCHKLL/MI//-.5 $ DIAGONAL TERM AND ROOT 

DIAGONAL KBAR/K1//.S $ DIAGONAL TERM AND ROOT 

ADD MI,KI/Q/.15915S//1 S CALCULATE FREQUENCIES 

MATPRN OCHKLL,KBAR,QL,Q,// $ 

$ 

PARAML CMPHIXZ//’TRAILER7l/S,N,NOZSET S NUM RITZ = NUM EIGENVECTORS 

S 

LABEL ENDBF S 

$ 

$ TRANSFORMATION MATRIX CREATION 

$ - BASED ON RUBIN-MACNEAL FOR FREE MODES (DOESN’T USE INREL) 

$ 

ALTER 807 S 

COND TRANSF.BFLEX S GO TO TRANSFORMATION DMAP 

JUMP ETRANSF $ JUMP AROUND IF NOT REQ. 

$ 

ALTER 808 S 

LABEL TRANSF $ AVOID INREL MODULE 

$ 

$ FREE BOUNDARY TRANSFORMATION 

$ 

IF(NORC>0) THEN S EXECUTE ONLY IF RC SET 

S 

$ CALCULATE SYSTEM FLEXIBILITY WHEN RITZ VECTOR HAVE BEEN 
S ORTHOGONALIZED WRT THE MASS MATRIX (EQN. (15) ABDALLAH) 

S 

SMPYAD QL,CMKXX,QL,,,/GKl/S////l////0 S INTER CALC 
SOLVE GKl,/GKlINV/3 S INVERT INTER CALC 

TRNSP QL/QLT $ TRANSPOSE qL 

SMPYAD QL,GKlINV,QLT,,,/GK/3////////6 S CALC MODAL FLEX 

$ 

$ CALCULATE RESIDUAL FLEXIBILITY MATRIX 
$ NOTE: IF NO R-SET IS PRESENT THEN GE IS KINV 

S 

ADD GE,GK/GD//-1.0/0/ $ GD = GE - OK 

* 

% PARTITION DOWN SQUARE FLEXIBILITY MATRIX TO C OR BOUNDARY DOFS 

$ 

PARTN GD,VVRCO,/GCRCD,GORCD,,/-l//6/ S GD = (GCRCD/GORCD) 

$ 

$ FORM FREE * CONSTRAINT* MODE MATRIX, WHICH IS 
$ THE LOWER LEFT PARTITION OF THE TRANSFORMATION MATRIX 

S 

SOLVE GCRCD,/GCCDINV/3 $ INVERT GCRCD 

MPYAD GORCD.GCCDINV./GOATFRS/ $ GOAT = GOCD*GCCD**-l 

$ 

S FORM FREE * MODAL” MATRIX, WHICH IS 

S THE LOWER RIGHT PARTITION OF THE TRANSFORMATION MATRIX 
S 

PARTN CMPHIXZ„WRCO/QRCL,QOL„/l/ S CMPHIXZ = (QRCL/QOL) 

MPYAD GOATFRS,QRCL,QOL/PHIVZ//-l/ S PHIVZ=QOL-GOAT*QRCL 

$ 

S PACK GOATFRS WITH ZEROS TO MAKE IT A-SIZE 

S THE LOWER LEFT PARTITION OF THE TRANSFORMATION MATRIX 
$ 

VEC USET/VAQRC/’AV’Q’/’COMPy S FORM PART VECTOR 
MERGE, ,,GOATFRS )J VAQRC,/GOATFREE/I/ S MAKE A-SIZE 

S 

S FORM REPLACEMENT BOUNDARY MASS AND STIFFNESS MATRIX, THIS 
S TRANSFORMATION IS THE SIMILAR TO STANDARD CMS AND SOL 63, 

$ EXCEPT FOR FREE, REVISED ’GOAT’ MATRIXIS NOT EQUAL TO 
$ CONSTRAINT MODES, SO STIFFNESS TRNASFORMATION IS SAME AS MASS 

$ 

MPYAD MOO,GOATFREE,MOA /MOA1FR $ IDENTICAL TO 820 

MPYAD MOA,GOATFREE,MAAl/MA A2FR/I $ INTERMEDIATE 

MPYAD GOATFREE,MOAIFR,MAA2FR/MA A/1////0 S NEW MASS MATRIX 

$ 

MPYAD KOO,GOATFREE,KOA/KOAlFR S STIFFNESS SECTION 

MPYAD KOA,GOATFREE,KAAl/KAA2FR/l S INTERMEDIATE 

MPYAD GOATFREE,KOAIFR,KAA2FR/KAA/l////6 * NEW KAA MATRIX 
S 

$ REPLACE PREVIOUS GOAT WITH FREE ONE CALCULATED HERE 
$ 

EQUIV GOATFREE, GO AT/ALWAYS % 

S 

ENDIF $ 

S 

$ MAKE PHIVZ Q-SIZE, THIS REPLACES INREL FOR FREE AND FIXED 
S BOUNDARY FLEXIBILITY METHODS 
$ 

IF (NOZSET > NOQSET) MESSAGE / 
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/’FATAL ERROR-NUMBER OF MODES EXCEEDS Q-SET DOF’/f 
IF (NOZSET > NOQSET) JUMP RFERR $ 

$ 

MATGEN ,/ZQ/«/NOQ9ET/NOZSET/NOQSET S CREATE PARTITIONING MATRIX 
MERGE PHIVZ„„ZQ,/GOQ/l $ TO MAKE PHIVZ Q-SIZE (GOQ) 

$ 

COND PRINT#,DIAG S 

MATPRN CMPHIXZ,GKl,GKiINV > GK l GD//S 
MATPRN GCRCD,GORCD,GCCDINV,GOATFRS f QRCL//S 
MATPRN QOL,PHIVZ,GOATFREE,MAA,KAA// S 
LABEL PRINT# $ 

$ 

LABEL ETRANSF $ 

$ 

$ SETUP INTERMEDIATE MODAL STIFFNESS MATRIX 

S 

ALTER 616,616 S 

UMERGEI USET i KQQ m ,/KLAA1/ , AV»Q7*T70 $ 

$ 

$ CREATE OFF-DIAGONAL MODAL STIFFNESS TERMS 
S 

ALTER 816 S 

IF(NORC>0 AND BFLEXCO) THEN S 
MPYAD KOO,GOAT,KOA/KOAl S 
MPYAD GOAQ,KOAl,/KQT/l $ 

TRNSP KQT/KTQ $ 

ADD KQT,KTQ/KAAQT $ 

ADD KLA A1,KAAQT/KLA A $ 

ELSE $ 

EQUIV KLA Al, KLA A/ ALWAYS $ 

ENDIF $ 

$ 

ALTER 669$ 

COND PRINTT,DIAG S 
MATPRN GOQ,KLAA,MLAA„// $ 

LABEL PRINT7I 
S 

END ALTER 8 
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